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Abstract 

It has been recently shown that excess events observed by the LSND and Mini- 
i-C \ BooNE neutrino experiments could be interpreted as a signal from the radiative 

decay of a heavy sterile neutrino produced in neutral current-like neutrino in- 
^■i teractions. If the Vh exist, it would be also produced by the beam from the CERN 

SPS in the neutrino beam line shielding. The v^s would penetrate the shielding and 
be observed through the Vh —> 7^ decay followed by the photon conversion into e + e~ 
pair in the active target of the NOMAD detector. The Uh's could be also produced 
in the iron of the magnetic spectrometer of the CHORUS detector, located just in 
front of NOMAD. Considering these two sources of VhS we set new constraints on 
Vh properties and exclude part of the LSND/MiniBooNE parameter space using 
bounds on single photons production in neutrino reactions recently reported by the 
NOMAD collaboration. We find that broad bands in the parameter space are still 
open for more sensitive searches for the i>h in future neutrino experiments. 
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Over the past 10 years there is a puzzle of the 3.8 a event excess observed by 
the LSND collaboration [1] . This excess originally interpreted as a signal from 
Wfj, V e oscillations was not confirmed by further measurements from the 
similar KARMEN experiment [2]. The MiniBooNE experiment, designed to 
examine the LSND effect, also found no evidence for — » v e oscillations. How- 
ever, an anomalous excess of low energy electron-like events in quasi-elastic 
neutrino events over the expected standard neutrino interactions has been ob- 
served [3]. Recently, MiniBooNE has reported new results from a search for 
TJfj, — > V e oscillations |4j. An excess of events was observed which has a small 
probability to be identified as the background-only events. The data are found 
to be consistent with — > V e oscillations in the 0.1 eV 2 range and with the 
evidence for antineutrino oscillations from the LSND. 
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In the recent work [5] (see also [6fTf8] ) it has been shown that these puzzling 
results could all be explained in a consistent way by assuming the existence 
of a heavy sterile neutrinos (v^). The is created in neutral- current (NC) 
interactions and decay subsequently into a photon and a lighter neutrino v 
in the LSND and MiniBooNE detectors, but it cannot be produced in the 
KARMEN experiment due to the high energy threshold. The could be 
Dirac or Majorana type. The Uh could decay dominantly into 71/ pair if, for 
example, there is a large enough transition magnetic moment between the 
and v mass states. Assuming the is produced through mixing with u^, the 
combined analysis of the LSND and MiniBooNe excess events suggests that 
the Uh mass, mixing strength, and lifetime are, respectively, in the range 

40 < m h < 80 MeV, lO" 3 < \U, h \ 2 < lO" 2 , lO" 11 < r h < lO' 9 s. (1) 

A detailed discussion of consistency of these values with the constraints from 
previous searches for heavy neutrinos [9] is presented in [5] . Briefly, the mixing 
of ([1]) is not constrained by the limits from the the most sensitive experiments 
searched for extra peaks in two-body 7r, K decays PfTU] . because the v h mass 
range of ([TJ is outside of the kinematical limits for 71^2 decays, and not accessi- 
ble to K^2 experiments due to experimental resolutions. The parameter space 
of ([1]) cannot be ruled out by the results of high energy neutrino experiments, 
such as NuTeV [TlLj or CHARM [T2"iri3] . as they searched for z/^'s of higher 
masses {m h > 200 MeV) decaying into muonic final states (/ivri/, /i/zi/, fieu, ..) 
[9], which are not allowed in our case. The best limits on (C/^j 2 derived for 
the mass range ([T]) from the search for — > e + e~v decays in the PS191 ex- 
periment [H], as well as the LEP bounds [15], are found to be compatible 
with ([I]) assuming the dominance of the Vh decay. New limits on mixing (t/^l 2 
obtained by using the recent results on precision measurements of the muon 
Michel parameters by the TWIST experiment [TB] are also found to be con- 
sistent with (pQ). Finally, the most stringent bounds on (f/^hl 2 coming from the 
primordial nucleosynthesis and SN1987A considerations, as well as the limits 
from the atmospheric neutrino experiments, are also evaded due to the short 
Vh lifetime. 

Recently, several constraints on the properties of the obtained from the 
muon capture on hydrogen [XT] and radiative decays of charged kaons [TBTTQ"] 
have been reported. The new bounds, found to be in tension with ([1]), are 
obtained under assumption that the z/^, as a component of the z/ M , could 
be also produced in reactions induced by charge-current (CC) interac- 
tions. However, one may argue that the heavy neutrino explanation of the 
LSND/MiniBooNE anomaly assumes that excess events are originated from 
the decays of zVs produced in a new i/^NC-like interactions, i.e. without charge 
muon emission [5]. Therefore, the constraints of Ref.'s [TT|18|ll9j can be evaded, 
if e.g. the new fermion has new dominant NC-like interactions, while its 
CC reactions are suppressed. 
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Fig. 1. Schematic illustration of the experiment to search for the v^, — > decay 
of a sterile neutrino. If the v\ t 's exist, they would be produced in v„ interactions in 
the CERN WANF v„ neutrino beam line shielding. The particles would penetrate 
the downstream shielding and be observed in the NOMAD neutrino detector through 
their — > decays into a photon and a light neutrino followed by the photon 
conversion into e + e~ pair in the NOMAD target. The shaded region of the shielding 
serves as a dump of the thickness about 9 \i n t ( interaction lengths ) to absorb neu- 
trino shower components which could escape absorption in the CHORUS detector 
matter, and generate background events in NOMAD. 

There are several possibilities for such scenario. For example, the can be 
produced by the Primakoff conversion, + Z — > + Z of the muon neu- 
trino in the electromagnetic field of nuclei due to dipole transition moment 
between the v h and states, with the subsequent decay u h — > due to 
the same interaction, see, e.g. Fig. 1 in |20j and, also [21 J for more detail 
discussions. Another mechanisms, we are mostly interested in, adopt new in- 
teraction to produce sterile neutrinos in muon neutrino scattering off nuclei. 
In the standard model, in our case of low energies, both Z- and VF-boson 
contributions to neutrino scattering are naturally of the same order, because 
there is no hierarchy between Z- and VF-boson masses. Would P^-boson be 
much heavier, charge current contributions to neutrino scattering off nuclei are 
suppressed as compared to that of neutral current. Consequently, to diminish 
the meaning of charge current channel in searches for sterile neutrinos one 
can involve a similar hierarchy for the new interaction responsible for sterile 
neutrino production. Such a model based on adding to the SM one new scalar 
doublet with the same gauge quantum numbers as the SM Higgs field [22] is 
discussed elsewhere [21]. Similar scenario could arise also in the low-energy 
superstring-inspired Eq theories [23]. In all these cases one does not need any 
sterile-active neutrino mixing at all to explain sterile neutrino appearance in 
muon neutrino beam traveling through a neutrino detector volume. Therefore, 
it would be interesting to obtain direct experimental constraints on the such 
v h properties. 

In this Letter we study new constraints on properties of iVs from the mea- 
surements of the NOMAD experiment at the CERN WANF neutrino beam, 
assuming the i/^s are produced dominantly in new z/^NC-like reactions. The 
present analysis as well as the experimental signature of the signal events are 
similar to those of heavy neutrino searches previously reported in [2^1l25|26] . 
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The CERN West Area Neutrino Facility (WANF) beam line [27] provides an 
essentially pure v^ beam for neutrino experiments. It consists of a beryllium 
target irradiated by 450 GeV protons from the CERN SPS. The secondary 
hadrons of a given sign are focused with two magnetic elements located in 
front of a 290 m long evacuated decay tunnel. Protons that do not interact 
in the target, secondary hadrons and muons that do not decay are absorbed 
by a 400 m thick shielding made of iron and earth. The NOMAD detector 
is located at 835 m from the neutrino target. The detector is described in 
detail in Ref . [2B] • It consists of a number of sub-detectors most of which are 
located inside a 0.4 T dipole magnet with a volume of 7.5x3.5x3.5 m 3 : an 
active target of drift chambers (DC) with a mass of 2.7 tons (mainly carbon), 
and a total thickness of about one radiation length followed by a transition 
radiation detector, a preshower detector, and an electromagnetic calorimeter. 
A hadron calorimeter and two muon stations are located just after the magnet 
coils. A plane of scintillation counters, V, in front of the magnet was used to 
veto upstream neutrino interactions and muons incident on the detector. 

If the Vh exists, it would be produced in interactions 

Up + N^Vh + X. (2) 

of muon neutrinos in the WANF neutrino beam shielding. If the v^, is a rela- 
tively long-lived particle, the flux of v^s would penetrate the downstream part 
of the shielding without significant attenuation and be observed in NOMAD 
through Vh —> 7^ decays followed by the decay photons conversion into e + e~ 
pairs in the NOMAD DC target, as schematically illustrated in FigJTJ The ex- 
perimental signature of v h — > jv decays would appear as an excess of isolated 
e + e~ pairs above those expected from standard neutrino interactions. As the 
NOMAD searched for an excess of isolated e + e~ pairs in their detector, the 
results obtained can be also used to constrain the radiative decay of the Vh- 
Note, that the search for an excess of isolated photon events in neutrino inter- 
actions is complicated by the presence of a large background coming from the 
electromagnetic decays of neutral mesons (mostly 7T°,t],t]' — > 77). However, 
because the experimental signature of these events is clean, they can be se- 
lected with significantly suppressed background due to the excellent NOMAD 
capability for precise measurements of e + e~ pairs, see for example [24"|25f26] . 

To make quantitative estimates, we performed simplified simulations of the 
Vh production (J2J) and decay in the NOMAD detector schematically shown in 
Fig. [TJ The flux of v^'s produced in the shielding was calculated by using the 
known flux of v^ beam in NOMAD [29J . In these simulation the Vh production 
vertex was required to be located in the upstream part of the shielding exclud- 
ing shaded region shown in Fig{TJ This end cap region serves as an additional 
dump to absorb neutrino hadronic shower components, which might escape 
absorption in the matter of the CHORUS detector located downstream. The 
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Fig. 2. Energy spectra of muon neutrinos that produce heavy neutrino with the mass 
of 50 MeV in neutral current interactions in the WANF neutrino beam line shield- 
ing. 

choice of the dump thickness of 9 Xi nt (interaction lengths) allows to fully 
protect NOMAD from the leak of secondary particles produced in the up- 
stream neutrino interactions, and suppress background from punch-through 
K\ and n interactions in NOMAD that could produce fake signal. More ac- 
curate analysis would require detailed simulations of the background in the 
NOMAD detector, which is beyond the scope of this work. Note, however, 
that obtained limits are not very sensitive to the choice of the dump thickness 
within the range ~ 9 — 15 X int . 

The simulated energy spectrum of the incident v^s is shown in FigfJJ It is 
assumed that, similar to the case of neutrino mixing, the cross section of 
heavy neutrinos production fl2]) can be expressed in a model-independent way 
as follows: 

a(v„ + N^u h + X) = a^oiy^ + N -» u„ + X)f ph , s , (3) 

where aiy^ + N — > + X) is the cross section for z/ M NC interactions, f p h. s . = 
a(m)/a(0) is the phase space factor which takes into account dependence on 
the v h mass, and a^h play a role of an effective coupling strength in new ivNC- 
like interactions. The distribution of energies for u^'s with momenta pointing 
to the NOMAD fiducial area, and photons from the isotropic — > decay in 
the NOMAD target are shown in Figj3j Once the flux was known, the next 
step was to calculate the e + e~ spectrum based on the v h — > •yu decay rate. 
For a given flux ^(uh), the expected number of signal events from — > yv 
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decays occurring within the fiducial length L of the NOMAD detector located 
at a distance V from the production vertex is given by 



N Uh _^ u = [ A<S>(v h )exp(-^^)[l-exp(-^)}^P,e e+e ~dE Uh dV (4) 

Vvh^h Vvy^h tot 

where p Uh is the v% momentum, is its lifetime at the rest frame, T e + e -, T tot 
are the partial and total z/^-decay widths, respectively, P 7 is the decay photon 
conversion probability, and e e + e - is the e + e~ pair reconstruction efficiency The 
acceptance A of the NOMAD detector was calculated tracing zVs produced in 
the shielding to the detector. As an example for a mass m Uh = 50 MeV, A = 
4.8% and e e + e - ~ 25%. It is assumed that the total rate Y tot of the v h decays 
is dominanted by the radiative decay v h — > 71/, see discussion in [5] , hence the 
branching fraction of the — > 71/ decay is BR{vh — > 71/) = r ^p~ 5>7 ^ ~ 1. 




Fig. 3. Energy distributions of Vh neutrinos with the mass 50 MeV produced in the 
v^NC interactions in the beam dump and pointing towards the NOMAD fiducial 
area, and photons from the — > 7^ decay in the NOMAD target. The spectra are 
normalized to a common maximum. 

The NOMAD search for an excess of single photon events is described in detail 
in Ref. [24]. Briefly, it used data collected with of 5.1 x 10 19 protons on target. 
The strategy of the analysis was to identify e + e~ candidates by reconstructing 
in the DC isolated e + e~ pairs with the total energy greater than 1.5 GeV and 
invariant mass below 100 MeV, that are accompanied by no other activity 
in the detector. The measured rate of e + e~ pairs was then compared to that 
expected from known sources. Two samples of signal photon candidate events 
were selected in NOMAD among the total number of 1.44 xlO 6 i^CC reaction 
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Fig. 4. Distribution of variable E^(l — cos0^ 7 ) for photons from the decay — > 
of the Vh with the mass 50 MeV produced in the u^NC interactions in the SPS 
neutrino beam dump and decaying in the NOMAD fiducial volume. 

recorded [21]. The first sample corresponds to photons emitted at any angle 
0jy 7 with respect to the primary neutrino beam direction, while the second one, 
corresponding to photons produced in the forward direction, was selected by 
using the cut ( < 0.05 GeV on variable ( = E^(l — cosQ ul ). After applying all 
selection criteria, 155 and 78 candidate events, with a predicted background of 
129.2±8.5±3.3 and 76. 6± 4.9 ± 1.9, were observed in two regions, respectively. 
These results are found to be consistent with the background expectations and 
yield an excess of 25.8 ±15.5 and 1.4 ±10.3 events, respectively. The measured 
spectra are also found to be in good agreement with prediction. Hence, no 
evidence for an excess of single photons produced in neutrino interactions 
has been observed and the corresponding upper limits of < 51 events and < 
18 events at 90% CL for the number of single photon excess events in each 
sample were obtained, respectively. 

Using Eq.fll]) and these upper limits, we can then determine the 90% CL 
upper limit for the corresponding coupling a^h- The distribution of variable ( 
for decay photons in NOMAD is shown in FigJH As the zVs arrive at NOMAD 
from the source located at a relatively far average distance of ~ 70 m, most of 
the decay photons are produced in the forward direction, at a small values of 
( and angle ©„ 7 . In this case, the best bounds are obtained by selecting signal 
events with the cut ( < 0.05 GeV. The corresponding 90% CL exclusion region 
in the a^h vs T h plane shown in Fig. [5] together with the result obtained with 
the CHORUS target (see below) is calculated by using the relation N Uh _^ lu < 
18 events. Our result is sensitive to a coupling a^h ^ 10~ 3 , and corresponds 
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Fig. 5. The 90% CL exclusion regions in the (t^ol^)- plane obtained from the 
NOMAD results [2^] for the case of the production in the SPS neutrino beam 
shielding and in the iron of the muon spectrometer of the CHORUS detector (Fe 
target). The area of the allowed LSND/MinibooNE parameter space corresponding 
to the Uh mass of 50 Me V is also shown. 

to the Vh lifetime region Th ~ 10~ 10 — 10~ 9 s. Over most of this region, the Vh 
lifetime is sufficiently long, that Lnih/phTh <C L'mh/VhJh "C 1. 

For the lifetimes shorter than ~ 10~ 9 s, new constraints on properties 
can be obtained by considering the Vh production in the muon spectrometer 
matter (mainly iron) of the CHORUS detector located just about 10 m up- 
stream of the NOMAD at the same neutrino beam line, as shown in FigJTJ 
The CHORUS detector, specifically designed to search for — v T oscillations 
at the SPS neutrino beam, is a hybrid setup which combines a nuclear 
emulsion target with various electronic detectors such as trigger hodoscopes, 
a scintillating fibre tracker system, a hadron spectrometer, electromagnetic 
and hadronic calorimeters, and the muon spectrometer [3U]. The muon spec- 
trometer consists of seven tracking sections interleaved by six magnets. Each 
magnet is made from twenty iron disks, each of 3.75 m in diameter and 2.5 cm 
thick, with scintillator planes interspersed. The total weight of this Fe-target 
is about 260 t. 



Similar to above considerations, the expected number of heavy neutrino is h — > 
7Z/ decays in the fiducial volume of the NOMAD detector can be estimated 
by using expression (j3J), and taking into account the corresponding number 
of heavy neutrinos produced in the primary neutral current interaction in 



8 



the CHORUS iron target. In this estimate we required the Uh production to 
occur in the first three upstream sections of the CHORUS muon spectrom- 
eter, assuming that the hadronic showers from the neutrino interaction are 
completely absorbed in the remaining three downstream sections served as a 
"beam dump" with the total thickness of 9 \ n t. This requirement is necessary 
to suppress background from K\ and n interactions discussed above, and also 
to avoid rejection of signal event by the NOMAD VETO system, when a sec- 
ondary energetic hadron trigger the NOMAD veto counter V. The simulations 
show that due to the shorter distance to the Fe target the angular distribution 
of the decay photons converted in the NOMAD DC target is wider in com- 
parison with the case, considered above, and the most restrictive limits are 
derived by using the NOMAD 90% CL upper limit of < 51 events obtained 
for the case of the wider ( distribution of the single photon events. 

The final 90% CL upper limit curves in the (r^; a^h) plane are shown in Fig. 
[5] together with the LSND-MiniBooNE Vh parameter space calculated for the 
Vh mass of 50 MeV. The attenuation of the Vh flux due to absorption in the 
shielding is found to be negligible. For example, our results can also be used 
to constrain the magnitude of the transition magnetic moment (fih) of the 
hypothetical u h , see for more details Ref. [5]. The estimated attenuation of the 
Uh - flux due to Vh electromagnetic interactions with matter in this case is found 
to be negligible, since for the Vh lifetime in the range 10~ n < t/j < 10~ 9 s, the 
fih ranges from ~ 10 -8 /ib to ~ 10 _9 /Xb (here \xb is the Bohr magneton), and 
the heavy neutrino mean free path in iron is more than 10 km, as compared 
with the iron and earth shielding total length of 0.4 km used in the SPS 
neutrino beam. 

In summary, in this work we study possible manifestations of the radiative 
decays — > 71/ of a heavy sterile neutrino in the CERN WANF neutrino 
beam. The existence of the was proposed for the explanation of the anoma- 
lous excess events observed in the LSND and MiniBooNE experiments. It is 
assumed that u^s are dominantly produced in new i/^NC-like reactions in or- 
der to evade recent constraints from the Uh searches in reactions induced by 
z/^CC interactions. We study the v h production occurring either in the CERN 
WANF neutrino beam line shielding, or in the iron of the muon spectrometer 
of the CHORUS detector located upstream of the NOMAD and followed by 
the radiative decay Vh l v with the subsequent conversion of decay photons 
into e + e~ pairs in the NOMAD active target. Using sensitive limits from the 
NOMAD experiment on single photon production in neutrino interactions, we 
derive new constraints on the properties. The obtained results allow to ex- 
clude zVs of ([I]) for the lifetime values 3 x 10~ n < r h < 5 x 10~ 10 s, while 
still leaving open two bands for the existence of Vh with the lifetime either 
T/i ~ 10 -9 s or Th < 3 x 10~ n s. We have demonstrated a significant potential 
of the discussed experimental approach for the future more sensitive searches 
for the Vh at the high intensity neutrino facilities. For example, more restric- 
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tive constraints are expected to be obtained from the search for the — > 
decays from the produced in z/ M interactions in the future LAr experiment 
[2Tj . aiming to check the origin of the LSND/MiniBooNE excess events at a 
neutrino facility at CERN [5T1I52] . 
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M.M. Kirsanov, N.V. Krasnikov, L. Di Leila, V.A. Matveev, F. Pietropaolo, A. 
Rubbia, C. Rubbia, and S. Mishra for their interest, and valuable discussions 
and/or comments. The help of D.S. Gorbunov and M.M. Kirsanov is greatly 
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